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ABSTRACT 

We investigate the physical properties of tidal structures in a disk galaxy 
created by gravitational interactions with a companion using numerical A-body 
simulations. We consider a simple galaxy model consisting of a rigid halo/bulge 
and an infinitesimally-thin stellar disk with Toomre parameter Q « 2. A per- 
turbing companion is treated as a point mass moving on a prograde parabolic 
orbit, with varying mass and pericenter distance. Tidal interactions produce 
well-defined spiral arms and extended tidal features such as bridge and tail that 
are all transient, but distinct in nature. In the extended disks, strong tidal force 
is able to lock the perturbed epicycle phases of the near-side particles to the 
perturber, shaping them into a tidal bridge that corotates with the perturber. 
A tidal tail develops at the opposite side as strongly-perturbed, near-side parti- 
cles overtake mildly-perturbed, far-side particles. The tail is essentially a narrow 
material arm with a roughly logarithmic shape, dissolving with time because of 
large velocity dispersions. Inside the disks where tidal force is relatively weak, 
on the other hand, a two-armed logarithmic spiral pattern emerges due to the 
kinematic alignment of perturbed particle orbits. While self-gravity makes the 
spiral arms a bit stronger, the arms never become fully self-gravitating, wind up 
progressively with time, and decay after the peak almost exponentially in a time 
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scale of ~ 1 Gyr. The arm pattern speed varying with both radius and time con- 
verges to Q — k/2 at late time, suggesting that the pattern speed of tidally- driven 
arms may depend on radius in real galaxies. Here, Q and k denote the angular 
and epicycle frequencies, respectively We present the parametric dependences of 
various properties of tidal features on the tidal strength, and discuss our findings 
in application to tidal spiral arms in grand-design spiral galaxies. 

Subject headings: galaxies: spiral — galaxies: structure — galaxies: interactions 
— galaxies: evolution — methods: numerical 



1. Introduction 



Spiral arms are the most outstanding morphological features in disk galaxies. They not 
only provide information on the dynamical states of the background stellar disks but also af- 
fect galactic evo l ution by triggering large-sc ale star formation in the gaseous component (see 



Elmegreenlll995l : lBertin fc Linlll996l ; see also lMcKee fc Ostrikerll2007l and references therein). 



Regarding the nature of the spiral structure, two pictures have been proposed. In one pic- 
ture, the arms are viewed as quasi- stationary density wav es that live long, rotating almost 



rigidly around the galactic centers (ILin &: Shulll964j . Il966l ). Nonaxisymmetric instability of 



the st ellar disks may grow to form a sort of se lf-sustained standing density waves in stellar 
disks faertin et al.lll989al lbl: [ Bertin &: Linlll996l ). In the other picture, the arms are t ransient 



1969 




1966 





Toomre fc Toomrelll972h or by s wing amplifica tion of leading waves ( Julian fc Toomre 



Goldreich fc Lynden-Belllll965l : iToomrd 119811 ). In this case, spiral features are short 



lived, lasting only for several rotation periods (~ 1 Gyr) and perhaps requiring intermittent 



external forcing (e.g.. ISellwood &: Carlbergill984j ). 



Observations indicate that the probabilit y to have grand- design arms is much higher for 

galaxies in binaries or groups than in the field ( Kormendy fc NormanlE979 : Elme green fc Elme green 
19821 . 119871 ). This suggests that regardless of their nature, some of grand-design spiral arms 
are clearly excited by nearby galaxies through tidal interactions. Prototypical examples in- 
clude M51 and M81 that possess, respectivel y, companion galaxi e s NG C5195 and M82 within 
50 kpc in distance. In a pioneering work, iToomre fc Toomrei (Il972l ) used non-interacting 
test-particle simulations to demonstrate that gravitational interaction of a disk galaxy with 
its companion generates features such as tidal bridge and tail which are commonly seen in 
extended disks of inter acting galaxies. Inclusion of self-gravity tends to enhance spiral struc- 



ture in the disks (e.g., iHernquistl 1 1 990l ) . Grand-design spiral arms can be produced even by 



a low-mass perturber if the interaction involves a very close passage, indicating that tidal 
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arms may be more frequent than previously thought (IByrd fc Howard! 1 1992T ) 



Since iToomre &: Toomrd (Il972l ). there have been many numerical studies of tidal in- 

ogies in the 



teractions of galaxies, incl uding the detailed modelings of the arm morpho 



M51/NGC5195 system (e.g. , lHernquistl[l99ol : Irloward & Byr3l99g 



2000aB burrell et al1l2003l ) and in the NGC 7753/7752 sy stem JSalo fc Laurikainen 



formation of tida^ 



Elme green et al 



tails and tidal dwarf galaxie s therein (e.g.. lBarnes fc Hernquistl 



Barnes 



1998 



Salo fc Laurikainen 



1992 



993), 



1996 



1993 



Wetzstein et al.l 120071). a nd formation of bars at the central parts of 



galaxies (e.g., Noguchi 1987 ; iGerin et aL 199oh . In particula r , by u sing numerical simula 
tions with self-gravitating stars and gas, ISalo &: Laurikainen! (j2000al ) argued that a bound 
multiple-passage orbit of NGC519 5 better reproduces the observed kinematics of an extended 
H I tail of M51 (IRots et al.lll990l ). whereas the radial velocity data of the pl anetary nebu- 



lae a ssociated with the tidal structures favor an unbound single-passage orbit (IDurrell et al. 



20031 ). Although these authors considered primarily the outer, extended tidal features for 
comparison, the spiral arms in the main disk may more tightly constrain the orbital pa- 
rameters of the M51 system. This is because the structure and kinematics of tidal tails 
depend rath er sensitively on the o bservationally uncertain parameters such as a halo mass 
distribution (IDubinski et al.lll996l ). 



While the aforementioned work has improved our understanding of the tidally-induced 
morphological changes of galaxies chiefly in extended disks, these studies did not focus on 
spiral structures in the main disks that are more relevant to large-scale star formation. 
Stellar spiral arms are certainly one of the main agents that greatly influence dynamical 
evolution of the interstellar gas in disk galaxies. Since the turbulent and thermal sound 
speeds of the gas are small, the gas responds very strongly to the gravitational potential per- 
turbations impose d by the stellar spiral arms, readily forming galactic spiral sho cks near the 
potential minima (jRoberts!ll969l ; IShu. Milione. fc Robertall973uWoodwardlll975l ). In optical 
images, these shocks appear as narro w dust lanes that represent regions where giant molec 



ring inside the dust lanes (e.g., 



ular clouds and new s t ars form fe.g.. lElmegreen fc Elmegreenl Il983l ; IVogel et all Il988l ; iRand 



19931 ; lElme green! 1 1 9 94 ; IShetty et al.ll2007l). Nonaxisymmetric gravitational instability occur- 



Balbuslll988l : iKim fc Ostrikerl 12002 . 120061 : IShettv fc Ostriker 



20061 ) is most likely responsible for observed arm substr uctures including gaseous spurs (or 



feathers) that jut perpendicularly from the ar ms (e.g., 
20041 ; ICalzetti et alJliooSlLa Viene et al.ll2006h . 



Scoville et al. 2001; Willner et al. 



The strength of spiral shocks and their susceptibility to gravitational instability are 
strongly affected by the physical properties of stellar arms such as amplitude, pitch angle, 
pattern speed, etc, yet it is quite difficult to characterize them observationally. While the 
arm pitch angle can be determined relatively straightforwardly if the inclination of a galaxy 
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is known, it is challenging to measure the pattern speed unambiguously. For instance, the 
Tremaine fc Weinberg! (11984) method that has b een applied to the CO data of several grand- 



design spiral galaxies (e.g., IZimmer et al.l 120041 ) assumes, among others, that the pattern 



speed is independent of radius and that the molecular gas satisfies the mass conservation 



equation. In the case of M51, however, a recent study by IShetty et al.l (120071 ) shows that 
observed density and velocity profiles across the disk do not obey the continuity equation in 
any frame rotating at a fixed angular speed. Also, the strength of spiral arms determined 



Rix & Rieke 


1993; 


Patsis et al. 


2001) 



(e.g., IKix fc Kiekel 119931 ; IFatsis et al.l 120011 ). Given these observational uncertainties, it is 
desirable to run numerical simulations to pin down the arm parameters and thus access the 
connection between stellar arms and large-scale star formation in gaseous arms. 

Motivated by these considerations, we in this paper use numerical iV-body simulations 
to explore in detail the properties of stellar spiral arms resulting from tidal interactions. 
Since the parameter space is large, we consider an idealized galaxy model in which an 
infinitesimally-thin, two-dimensional, exponential stellar disk is immersed in a combined 
potential due to rigid halo and bulge. Self-gravitating particles comprising the disk respond 
to a point-mass perturber that passes on a prograde parabolic orbit in the same plane as 
the disk rotation; we vary the mass and pericenter distance of the perturber to study the 
situations with various tidal strength. The particles are not allowed to move out of the disk 
plane, and the effect of gas is ignored. A fully self-consistent treatment of the problem, 
using active halo and bulge as well as three-dimensional disks consisting of both stars and 
gas, will be stud i ed in subsequent papers. Similar simulations have been carried out by 
Elme green et al.l (Il99ll ) who showed that a cold stellar disk (with zero velocity dispersion) 
turns into a transient ocular shape if tidal perturbations are strong. In this work, we instead 
consider a disk galaxy with realistic velocity dispersions. Our main objectives are to study 
the quantitative changes in the properties of spiral arms as the tidal strength varies, and 
also to clarify the development and physical nature of tidal features known as bridge and 
tail in extended disks. 

This paper is organized as follows. Section [2] describes the galaxy model and the orbital 
parameters of tidal interactions as well as the numerical method we use. In §31 we focus on 
the transient extended-disk structures produced by strong tidal perturbations and show that 
the tidal bridge and tail form by distinct mechanisms. In §U we measure the properties (pitch 
angle, strength, and pattern speed) of the spiral arms and present their temporal and radial 
variations. Finally, we summarize our results and discuss their astronomical implications in 

m 
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2. Model and Numerical Method 

2.1. The Model Galaxy 

In this paper we investigate the generation of tidal features and their properties in a 
disk galaxy via gravitational interactions with a point-mass perturber. The disk galaxy 
consists of three components: a spherical "dark" halo, a spherical bulge, and an exponential 
stellar disk; we do not consider a gaseous disk in the current work. The halo and bulge 
accounting for the inner linearly-rising part and the outer nearly-flat part of the rotation 
curve are represented by fixed gravitational potentials for simplicity. This will ignore the 
potential consequences on the disk through the tidal deformation of the halo and bulge E On 
the other hand, this inert halo and bulge enables a large number of particles for the stellar 
disk. In order to maximize the particle number for stars near the disk midplane, we impose 
a constraint that the disk remains infinitesimally thin during its whole evolution. 

Appendix [A] describes the specific model we employ for each component of the galaxy: 
a truncated logarithmic potential for the dark matter haloE a Plummer potential for the 
spherical bulge, and an exponential density profile for the disk. The total galaxy mass of 
M g = 3.24 x 10 11 M inside R = 25 kpc is dominated by the halo; the disk takes 16% of 
the total. We realize the infinitesimally-thin disk by distributing N = 514, 000 equal-mass 
particles on the disk plane and by assigning to them random velocities corresponding to 
the Toomre parameter of Q ~ 2. This value of Q fairly well represents the stellar disk 
in the solar neighborhood and is large enough to prevent spont aneous generation of spiral 



arms via swing amp lification in the absence of tidal forcing (e..g.. ISellwood fc Carlbergill984l ; 
Bertin et al.lll989bl ). Before applying tidal perturbations, we evolve the galaxy in isolation 
for two Gyrs to relax the phase space distribution into a global equilibrium. Appendix |B] 
presents the temporal evolution of an isolated disk and radial profiles of various quantities 
when an equilibrium is reached. We take the particle distribution at 1 Gyr and use it as an 
initial condition for tidal encounter experiments. This guarantees that morphological and 
structural changes of the disk occurring during interactions with the perturber are entirely 
due to tidal perturbations. 



1 Since the velocity dispersions of dark matter particles are usually much larger than those of disk 
stars, the impact o f the pe rturber to the disk through the live halo and bulge is small, as confirmed by 
Salo fc Laurikainenl |2000ri ). 



2 We have also run models without halo truncation and checked that the properties of tidal features inside 
25 kpc are almost indistinguishable from those under the truncation. 
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2.2. Perturber and Model Parameters 



As a perturbing companion, we consider a point particle with mass M p that moves on 
a parabolic orbit relative to the center of the galaxy in a prograde fashion. To study the 
excitation of spiral arms as cleanly as possible (i.e. without disk warping and bending waves) 
and to be consistent with the thin-disk approximation, its trajectory is confined to the same 
plane as the galactic disk. Assuming that the galaxy whose center lies at R = is spherical, 
the relative orbit (Rp, <f> p ) of the perturber in the polar coordinates is given parametrically 
by 



Rr. 



R 



9p3 
^-"'peri 



■ HTl V 1 ) 1 

1/2 



(x + £ 3 /3), 



X 



(2) 



G(M g + M p )\ 

where M g is th e total galaxy mass with in 25 kpc, x = tan(0 p /2), and -R per i is the pericenter 
distance (e.g., IPress fc Teukolskyl 119771 ). Note that t = (or x = 0) corresponds to the 
pericenter passage of the perturber. 

To explore tidal encounters with various strength, we consider nine self-gravitating mod- 
els that differ only in the mass and the pericenter distance of the perturber. We also run 
one non-self-gravitating model to study the effect of self-gravity on the arm properties. Ta- 
ble [H lists the parameters of each model and some simulation outcomes. Column (1) labels 
each run. Columns (2) and (3) give the perturber masss relative to the total galaxy mass 
and the pericenter distance, respectively. Column (4) lists the dimensionless tidal strength 
parameter defined by 



R„ 



R r 



AT 



(3) 

L g / \-" , peri/ \ J- J 

which measures the momentum imparted by the perturber to a disk particle at R g = 25 kpc 
relative to its original angular momentum (jElmegreen et al.lll99ll ). Here, AT is the time 
elapsed for the perturber to move over one radian near the pericenter relative to the galaxy 
center, and T = (R^/GMg) 1 ^ 2 is the time taken by stars at R = R g to rotate one radian 
about the galaxy center. Columns (5) and (6) give the fractions of the disk particles that are 
captured by the companion and those escaping from the whole system, respectively. Column 
(7) gives the time t ta i\ when the tidal tail becomes strongest, while columns (8) and (9) list 
the pitch angle i tSu n and surface density E ta ii of the tail at t = £ ta ii- Finally, column (10) gives 
the peak strength of the spiral arms. Model A2* is identical to model A2 except that the 
self-gravity of density perturbations in the disk is artificially taken to zero in the former. Note 
that the self-gravitational potential of the unperturbed axisymmetric disk, as represented by 
equation (1A4[) . is still included in model A2* to make the rotation curve intact. Models Al 
and C3 correspond to the strongest and weakest encounters, respectively. 
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In our presentation, the units of length and velocity are 1 kpc and 1 km s _1 , respectively, 
which give the characteristic time unit of t = 0.98 Gyr. All the simulations run from 
t/t = —1.0, corresponding to (R p , <f) p ) = (178.3 kpc, —45.6°) for our fiducial model A2, to 
t/t Q = 3.0. Seen from the above, the perturber passes through the pericenter (i? peri ,0) at 
t = in the counterclockwise direction which is the same sense as the disk rotation. 



2.3. Numerical Method 



To evolve disk particles in response to tidal perturba tions, we use t he GADGET code 
that is parallelized on a distributed-memory platform (jSpringel et al.l l200ll ). In GAD- 
GET, the evaluation of gravitational force uses the Barnes-Hut hierarchical tree algorithm 
( jBarnes fc Hutlll986l ) and assumes a spline-softened mass distribution of a point mass. Ex- 
cept at the beginning of the simulations, GADGET employs a new cell-opening criterion 
Ml 4 > «|a i d |r 6 , which produces, at a lower computational expense, force accuracy compa- 
rable to that obtained from the standard criterion r > 1/9. Here, a and 9 specify prescribed 
force error tolerances, M and I are the total mass and size of a cell, r is the distance of a 
particle to the center-of-mass of the cell, and a Q id is the gravitational acceleration on the 
particle computed at the previous timestep. For all the models presented in this paper, we 
adopt a = 0.02 and 9 = 0.8. 

For the gravitational softening, we take a soft ening; length of h = .4 kpc; the equivalent 
Plummer softening length is e = h/2.8 = 0.14 kpc (jSpringel et al.ll200ll ). The relaxation time 
associated with the forc e softening am ounts to ~ cr 3 e/(7rG 2 E m), where a = ^Joro$ and 
m is the particle mass (IRybickil Il97ll ). Since this time is longer than 10 Gyrs for R > 1 
kpc and N = 5 x 10 5 , tidal feature s that form in the stellar disk are not contaminated by 
particle noises and relaxation (e.g.. IWhitelll988l ). Particles are advanced by a second-order 



leapfrog scheme with fully adaptive and individual timesteps. All the simulations have been 
performed on an IBM p690 cluster using 16 processors, taking typically ~ 25 hours for a 
single run. 



2.4. Limitations of This Work 

In this work we employ highly idealized models of galaxies, perturbers, and their tidal 
encounters, and consider a limited range of tidal strength. This obviously introduces a few 
important caveats that should be noted from the outset: 



1. An infinitesimally-thin stellar disk imposed in the simulations neglects non-planar mo- 
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tions of stars such as in vertical oscillations and warps. It also overestimates self-gravity 
at the disk midplane. 

2. While a perturbing companion more likely has an extended density profile in real 
situations, we represent it as a point mass, which may overestimate the tidal force at 
the closest approach, possibly affecting the shape and structure of tidal bridge and tail 
that form in extended disks. 

3. Since we treat the galactic halo and bulge as being dynamically inactive, it is convenient 
to evolve the entire system in the coordinates centered on the center of the galaxy. This 
naturally ignores centrifugal and Coriolis forces that arise from the orbital motion of 
the galaxy relative to the center of mass of the whole system. The neglect of the 
indirect forces may spuriously suppress the growth of m = 1 spiral modes in the stellar 



disk, where m is the azimuthal wavenumber (e.g., lAdams et al.l Il989l ; lOstriker et al. 



19921 ). although it is unlikely to much affect m = 2 and higher order modes. 



4. By employing the prescribed parabolic orbit for a perturber, we neither consider the 
back reaction of the stellar disk to the perturber nor allow multiple encounters that 
would occur if the perturber is in a bound orbit. Furthermore, the prescribed orbit and 
the rigid halo and bulge do not allow us to capture the potential effects of dynamical 
friction and ensuing orbital decay of the interact ing galaxies, which may make the tidal 
tails longer and stronger (see e.g.. iBarnesI 119881 ). 



5. Limited to the cases with S £ 0.25, tidal tails created in our models are relatively weak 
and survive only for ~ 0.3 Gyrs (see §3.21) . The current weak- or moderate-encounter 
models preclude the possibility of prominent tails found in many interacting systems 
that live long (~ 1 Gyr or longer) and sometimes fragment into tidal dw arf galaxies 
(e.g.. lBarneslll988l . 1 19921 ; iBarnes fc Hernquistlll992l ; IWetzstein et al.ll2007l ). which may 
occur when tidal interactions are very strong. 



Given these constraints and limitations, we by no means attempt to reproduce tidal 
deformation of real galaxies. We instead focus on the formation mechanisms and physi- 
cal nature of tidally- driven disk structures, and compare the simulation results with the 
predictions of analytic theories, for which the simplifications made above are appropriate. 



3. Extended Tidal Features 



Using a restricted three-body technique, iToomre fc Toomrd (119721 ) demonstrated that 
tidal perturbations distort the extended portions of a disk to produce elongated and narrow 
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features, phenomenologically termed "bridge" and "tail" . The bridge is built at the near side 
of th e disk toward the perturber, while the tidal tail or "counterstream" forms at the far side 
(e.g.. iPfleidererl Il963l ) . Self-consistent numerical simulations including the disk self-gravity 



show that tidal perturb ations excite not only extended tidal streams but also spiral arm s 



in the main disks (e.g., iHernquistl Il990l ; iByrd fc Howard! Il992l ; ISalo fc Laurikainenl Il993l ). 



In this section, we focus on extended tidal features and distinguish between the physical 
mechanisms that form bridge and tail, some of which have previously been overlooked. 



3.1. Tidal Bridge 

To illustrate the dynamical responses of extended disks to a tidal perturber, we begin 
by presenting in detail the results from our fiducial model A2 with M p /M g = 0.4 and 
-Rperi = 35 kpc. Evolution of the other models are qualitatively similar. Figure [1] shows the 
morphological evolution of the stellar disk in model A2. The arrow and the associated number 
in each panel indicate the direction and distance (in kpc) to the perturber, respectively. Only 
20% of the particles are plotted to delineate tidal features from the disk. Figure [2] displays 
the perturbed surface densities of model A2 in the <fi — \ogR plane. At early time when 
the perturber is far away from the galaxy (t £ — 0.1), the tidal deformation of the disk is 
vanishingly small. As the perturber approaches the pericenter, the disk begins to undergo 
significant morphological changes, first forming a bridge (t ~ 0.0 — 0.1) at the outskirts of 
the disk close to the perturber and then a tail at the opposite side (t ~ 0.2). 

Tidal force imposed by the perturber excites the epicycle orbits of individual particles. 
In Appendix [UJ we use an impulse approximation to estimate the amplitudes 5R of the 
perturbed epicycle orbits in an averaged sense. Figure [3] plots as thin lines the resulting SR 
with differing M p based on the impulse approximation. The thick line is for the case with no 
tidal perturbations in which the radial oscillations of particles are purely due to the initial 
velocity dispersions. Also plotted as various symbols are the dispersions ((R — Rq) 2 ) 1 ^ 2 of 
the particle positions R at t — with respect to the initial locations Rq for models A2, B2, 
and C2. Here, the angular brackets ( ) denote an average over the particles in a given radial 
bin. Note that the numerical results are in good agreement with the corresponding analytic 
estimates over a wide range of radii. In regions of disks with R ^ 15 kpc, the deviation 
from the original epicycle orbits is quite small. It nevertheless enables well-defined spiral 
structure there, as we will discuss in §U In the extended disks, on the other hand, strong 
tidal perturbations severely affect the orbits of particles, causing them to traverse over large 
radial distances. 

Since the tidal force is asymmetric, particles at the near side to the perturber are 
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more easily pulled radially outward and will subsequently find themselves subject to even 
greater tidal force at larger Particles whose velocities exceed the escape velocity become 
unbound, and eith er are captured by the p erturber or escape from the combined galaxy- 



perturber system (IToomre &: Toomrd Il972l ). The fraction of the captured particles and 
the non-captured, freely escaping particles are given in columns (5) and (6) of Table [TJ 
respectively; these are fitted roughly with M cap /M d = 0.95S 2 ' 08 and M esc /M d = 0.Q7S 2 92 
for 0.04 S £ 0.3. In model C3 with S = 0.029, the tidal force is too weak to accelerate 
particles to the escape velocity. Although less than 7% of the total even in our strongest 
encounter model Al, the amount of mass stripped off by the tidal force depends fairly steeply 
on the tidal strength, and can be substantial for encounters with large S. 

While the bridge is a pathway through which mass transfer occurs, it also contains a 
significant amount of bound particles. Due to strong tidal force, the orbits of these bound 
particles are eventually arranged in such a manner that the maximum radial velocities always 
occur in the direction to the perturber while the perturber remains close to the pericenter. 
This is well illustrated in Figure H] which plots the azimuthal distributions of the particle 
velocities at R = 20 kpc in model A2 for t < 0.3. In each panel, the vertical dotted lines 
indicate the direction, <p p , to the perturber. Although the morphological change of the disk 
is almost absent when t = —0.1 (see Fig. [1]), the signature of the tidal interaction is already 
apparent in the azimuthal variations of the particle velocities. Note that in the bridge both 
Vr and \dv^/d(j)\ are maximized at <fi = <f) p (t) at the epochs shown in Figure HI That is, 
the phases of particle orbits are locked to the perturber during this time interval. Since the 
epicycle motions occur in the opposite sense to the disk rotation, this phase locking implies 
that v$ steadily decreases as the particles continue galactic rotation past the perturber. It 
attains minimum values near the leading edge of the bridge. It is at this leading edge where 
the particles fall rapidly radially inward, rendering the leading boundary of the bridge rather 
sharp. 

Figure [5] displays distortions of rings at several different initial radii Rq during the early 
phase of the tidal encounter. Near-side particles in a ring with larger Rq are pulled out earlier 
and by greater amount toward the perturber, shaping the ring into an egg-shaped oval. The 
tips of outer ovals become lagging behind the perturberQ At the same time, new particles 
from inner rings that rotate fast are pulled out to lead the perturber. This constructs a 



3 In the case of model A2, the ratio of the tidal forces at the near and far sides is 2.0 and 7.5 at R = 8 
and 20 kpc, respectively. 

4 The perturber in model A2 has an angular velocity of Q = 9.54 km s _1 kpc -1 at the pericenter. Since 
the corresponding corotation radius is R = 25 kpc in the disk, all the near-side particles shown in Figure [5] 
would lead the perturber were it not for strong tidal perturbations and the resulting phase locking. 
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transient pattern that persists while the perturber is close to the pericenter (~ a few tenths 
of Gyrs), with the pattern speed roughly equal to the instantaneous angular velocity of the 
perturber. As Figure [T] shows, the bridge in model A2 (and also in other models) lasts 
until t ~ 0.3 after which the perturber is too far away to tightly enforce the alignment of 
the epicycle orbits. Therefore, the bridge is a transient structure that not only allows mass 
transfer to the companion but also consists of bound particles that execute coherent forced 
oscillations in response to the applied tidal perturbations. 



3.2. Tidal Tail 

Tidal torque applied at the far side of the disk causes the leading (lagging) particles 
with respect to the line connecting the disk center and the perturber to lose (gain) angular 
momentum and thus to rotate slower (faster). This gives rise to a negative gradient of the 
circular velocity along the azimuthal direction. One may naively expect that the compressive 
velocity fields in the azimuthal direction should be a cause of a tidal tail at the far side, but 
this is not the case. The third panel of Figure H] shows that the velocity gradient amounts 
to dVffy/dcf) ~ —50 km s _1 rad^ 1 for the far-side particles at R = 20 kpc when the perturber 
is at the pericenter. Assuming that this value remains constant over a time interval At, the 
resulting fractional change <5£/£ of the surface density would be ~ Atdv^/^Rdcp) ~ 0.5 for 
At = 0.2 Gyr, which is too small to build a tidal tail in its own right. Therefore, the tail 
formation must involve additional processes. 

Figure [5] demonstrates the tail-making process in our models. Let us pay attention to the 
two groups of particles, denoted by dots in black or cyan, in the ring with R = 22 — 24 kpc. 
The dots in cyan representing a group of far-side particles at t = —0.05 are slowly rotating 
about the disk center, with a period of ~ 0.6 Gyr, by following moderately perturbed epicycle 
orbits. With relatively weak tidal force, the locking of the epicycle phases is not significant 
at the far side. On the other hand, the near-side particles in black that were ahead of the 
perturber at t — —0.05 have highly perturbed orbits, plunging toward the disk center as 
deep as R ~ 9 kpc at t = 0.05. The constraint of angular momentum conservation requires 
the particles to rotate faster at small R, providing them with a shortcut route to reach the 
far side (t ~ 0.10 — 0.15). A tidal tail develops as these strongly-perturbed, fast-rotating, 
near- side particles catch up with those mildly- perturbed, far-side particles (t ~ 0.15 — 0.2) 
(e.g. 



Pfleidererl 1 19631 ; iToomre fc Toomrd Il972l ) . 



Note that the outer-disk particles located in between the black and cyan dots, i.e., the 
particles with ~ 7r/4 — 7r at i = in the red ring in Figure [5], are all gathered into the tail 
extending to ~ 40 kpc from the disk center. Since the tail at a given radius is comprised of 
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particles from a wide range of radii in the unperturbed disk, it has large velocity dispersions 
both in the radial and azimuthal directions (t = 0.2 frame of Fig. H]). Accordingly, the tail 
in model A2 becomes weak and dispersed as the particles continue galactic rotation. This 
implies that the tails in our models are transient material arms. 

Figure [2] shows that the tail in model A2 forms at t ~ 0.2, is more or less logarithmic 
in shape with a pitch angle of tani ~ 0.5, and becomes more pronounced than the bridge. 
Time of tail formation £ tail , its pitch angle z ta ii, and its strength as defined by the surface 
density E ta o at R = 20 kpc and t = t tail of course depend on the the strength of tidal 
perturbations and the presence of self-gravity. Columns (7)-(9) in Tabled] list t ta ii, tanz ta ii, 
and E tai i for models with S > 0.07; models B3, CI, C2, and C3 with very weak tidal 
perturbations do not produce readily identifiable tail structure. These values are plotted in 
Figure E] as solid circles against S, which are well fitted with power laws: t tai \ = COTS' -0 54 , 
tani ta ii = O.TSS' ' 22 , and E t aii/S 20 = 79. OS 10 ' 72 , where S 20 indicates the surface density of the 
initial disk at R = 20 kpc. Definitely, a tail develops earlier and stronger for stronger tidal 
perturbations, although the pitch angle depends weakly on S. Note that the tidal tail in the 
non-self-gravitating model A2* is weaker and more loosely wound than that in model A2. 

As mentioned above, a tidal bridge in the near side consists of particles in coherent 
forced oscillations, while a tail in the opposite side forms by temporary particle overlapping. 
Both are transient features whose amplitudes decay after t ~ 0.2 — 0.3. As the perturber on 
a parabolic orbit moves away from the galaxy in our models, the diminished tidal force no 
longer aligns the phases of the particle orbits in the bridge. In addition, the large velocity 
dispersions of the tail are unable to keep it as narrow as when it first forms. Consequently, 
the particles making up the bridge and tail gradually spread out and follow galactic orbits 
with large eccentricities. They interact with each other and also with spiral arms, producing 
complicated structures seen at the extended parts of the disk in Figures [T] and [2J The further 
diffusion and interactions of particles eventually make the outer disk almost featureless in 
our simulations. 



4. Disk Structure 



We have seen in §3.11 that the enhancement of epicycle amplitudes due to tidal per- 
turbations is rather small in regions of disks with R £ 15 kpc. Nevertheless, the phases of 
perturbed epicycle orbits at different radii drift at different rates and are kinematically orga- 



nized to develop a trailing two- armed spiral pattern there (e.g.. lToomrdll969l ; iDonner et al. 
19911 ) . Figure [7] displays close-up views of density snapshots of model A2 in the x-y plane. 
A well-defined, two-armed spiral pattern is apparent for t ~ 0.2 — 1.0, becoming most con- 
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spicuous at t ~ 0.3 — 0.4. The spiral arms that appear as straight lines in the — lni? plane 
(see Fig. [2]) are approximately logarithmic, with a pitch angle varying with time. 

An inspection of Figure [2] reveals that the arms extend inward up to R ~ 4 kpc, 
corresponding to the inner Lindblad resonance (ILR), and smoothly join the the extended- 
disk features at t = 0.2. Since the pattern speed of the arms is smaller than the angular 
speed of the disk rotation, however, they soon decouple from the tidal tail (t = 0.3), and 
from the bridge at later time (t = 0.4 — 0.5) when the phase locking becomes inefficient. 
The spiral arms in our models are not stationary in the sense that their pattern speed is 
not constant over radius and that their pitch angle and amplitude vary with time. In this 
section, we explore the quantitative properties of the spiral arms. 



4.1. Pitch Angle 



Since the spiral arms that form in our models are logarithmic, it is useful to define the 
Fourier coefficients in 6 and In R as 



1 - 

A(m,p) = —^expiilmfa+plnRj}), 



(4) 



where N is the total number of particles, (Rj, (f)j) are the coordinates of the j-th par- 
ticle, and p is related to the pitch angle of an m— armed spiral through tanz = m/p 
( ISellwood fc Carlbergi Il984l ; ISellwood fc Athanassoulal Il986l ). A positive (negative) value 
of p corresponds to trailing (leading) spirals. 

Figure M plots the temporal evolution of the Fourier amplitudes |A(2,p)| of the m = 2 
logarithmic spiral mode in model A2 before the arms reach the maximum strength (t < 0.3). 
We consider particles only at R = 5 — 10 kpc where the pattern in model A2 achieves 
large amplitudes and contamination from the bridge and tail is almost absent. At early 
time (t £ 0.04), the modal growth occurs as the dominant p shifts from negative to positive 
values. This is suggestive of mild swing amplification in which seed perturbations grow as 
they change fro m leading to trailing. Since the corresponding amplification factor is less than 



10 when Q ~ 2 (IGoldreich fc Lynden-Belllll965l ; I Julian fc Toomrelll966 



Toomre 



1981) and 



since swing amplification becomes no longer efficient at p ^ 5 (e.g., ISellwood fc Carlberg 
19841 ). however, the further growth of the spiral modes cannot be attributable to swing 
amplification. It is rather due to the kinematic e ffects, enhanced by self-gravity, of the 
perturbed epicycle orbits in a manner described in iToomrd (119691 ). As the phases of the 
epicycle orbits drift and are coherently arranged, the density associated with the pattern 
grows quite rapidly and saturates at t ~ 0.3 in model A2. 
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It is well known that kinematic density waves without self-gravity tend to wind up due 
to the background differential rotation, with the pitch angle varying as 



tanz = t 



d(Q - k/2) 



dlnR 



(5) 



Binnev fc Tremainel 1987T). In the theory of quasi-stationary density waves hypoth- 



(e-g. 

esized by iLin fc Shul (11964 



19661 ). self-gravity of the spirals compensates for the winding 



tendency of the arms, keeping their pattern speed at a constant value over a wide range of 
radii. In order to check if this is the case in our simulations, we calculate the pitch angle of 
the arms determined from p that maximizes |A(2,p)| at a given time. Figure M shows the 
temporal changes in tani of the spiral arms located at R = 5 — 10 kpc, 8 — 13 kpc, and 11 — 16 
kpc for the self-gravitating models A2, B2, and C2, respectively. For comparison, Figure [9] 
also plots the results of the non-self-gravitating model A2* for the arm segments in an an- 
nulus with R = 8.0 — 8.5 kpc, over which d(Q — k/2) /dlnR is almost constant. Although 
the arm pitch angle in model A2* exhibits small fluctuations, the late-time portion can be 
well described by tanz oc consistent with the theoretical prediction (eq. [5])E For the 
self-gravitating models, the arms have moderate pitch angles amounting to tanz ~ 0.3 — 0.4 
when they grow and stand out initially. After attaining substantial strength, they begin 
to wind as tanz oc t~°- 5 °- 6 , with a smaller power index corresponding to stronger arms. 
This suggests that although self-gravity reduces the winding rate considerably, it cannot 
completely suppress the winding tendency of the spiral arms in our models^ 

Once finding the arm pitch angle and pattern speed (see below), we are able to compare 
the WKB theory of linear density waves with the simulation results. The local theory for 
tightly-wound linear waves in a stellar disk states that the perturbed radial velocity Sv r and 
azimuthal velocity 8vs are related to the perturbed surface density <5£ through 



Sv R = - — I — 



UK 

k R 



5E\ 



E J 



Svj, = - 



2Qk R F v {x) 



IT 



(6) 
(7) 



Since the Fourier method picks up, in a given annulus, the most dominant spiral modes that propagate 
radially inward, the time dependence of tan? can also be affected by the radial variation of d(fl — k/2) /din R 
if the annulus is wide enough. For instance, the average pitch angle of the arms in the R = 8 — 13 kpc region 
in model A2*, over which d(Cl — k/2) /dlnR varies by 13% relative to the mean value, decays as tan i oc i -0 94 . 

6 In addition to the background shear, short trailing waves in the presence of self-gravity would increase 
their radial wavenu mber kn as th ey propagate inward from the corotation radius, capable of decreasing the 



pitch angle further ( Toomre 1969h . 
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( ILin. Yuan, fc Shulll969l ). Here, = m/(Rtani) is the local radial wavenumber of the 
waves, x = (/c^Oft/ft) 2 , v = (Q p — mQ)/K is the dimensionless angular frequency with Q p 
denoting the patt ern speed, and F v {x) and are the reduction factors defined by equations 
(B9) and (B17) of ILin. Yuan. &: Shul (119691 ) . respectively. In equation ([7]), the imaginary unit 
i represents the phase shift between 5v^ and 5E. 

Figure [10] gives exemplary comparisons between the numerical results and the linear- 
theory predictions for the azimuthal variations of the perturbed variables. Two sets of 
numerical data near R — 10 kpc at t = 0.4 in model A2 and at t = 0.5 in model C2 are 
arbitrarily taken. In the top panels, black curves with some fluctuations draw from 
numerical simulations, while red lines plot the corresponding m = 2 Fourier modes <5£ m= 2. 
In the middle and bottom panels, red curves draw equations §6§ and ([7]) corresponding to 
5E m=2 . Blue curves represent the azimuthally-binned averages of v r and Sv^ = v $ — that 
are plotted as dots from the simulations. Apparently, the perturbed density in model C2 is in 
the linear regime and dominated by the m = 2 mode. Note that in spite of large dispersions 
in vr and Sv^, there is fairly good agreement between the numerical and analytic results 
for model C2. On the other hand, the spiral arms in model A2 are asymmetric and clearly 
in the nonlinear regime. I n this case, the per turbed velocities have significant contributions 



from high-m modes (e.g., IVandervoortl Il97ll ) . rising more steeply than a simple sinusoidal 



curve as particles leave the spiral arms. 

Among the models listed in Tabled] we found that models B3, C2, and C3 with relatively 
weak tidal perturbations (S < 0.06) produce linear spiral arms with sinusoidal density 
distributions. All the other models we considered show significantly nonlinear features in 
the density and velocity profiles. This implies that tidally-excited stellar spiral arms in 
grand-design spiral galaxies probably have non-linear amplitudes. 



4.2. Arm Strength 

One of the key parameters that directly influence gas flows in spiral galaxies is the 
strength of stellar spiral arms. Stronger spiral arms imply larger enhancement of gas density 
at the galactic shocks and hence more active star formation. To quantify the arm strength, 
we define 

_ 2nG5t m=2 rs v 

F = ' (8) 

where 5S m=2 denotes the amplitude of <5£ m=2 . Since the corresponding gravitational poten- 
tial perturbation is given by 5$ m = —2ixG5Y> m /{k R + m 2 R~ 2 ) 1 ^ 2 for a tightly- wound spiral 
in an infinitesimally-thin disk, F measures the gravitational force due to the spiral arms in 
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the direction perpendicular t o the arms relative to the the axisymmetric radial force Rfl 2 in 
the u nperturbed state (e.g. 
2002h . 



Roberts! 1196a IShu. Milione. fc Robertslll973l : iKim fc Ostriker 



Figure [TT] plots the radial variations of F for the arms averaged over the time interval 
At = 0.4 centered at the epoch when the arm amplitudes are maximized. In a given model, 
F broadly peaks at a certain range of radii; -R max ~ 5 — 10 kpc for models Al, A2, Bl, and 
CI, -R max ~ 8 — 13 kpc for models A2*, A3, and B2, and i? max ~ 11 — 16 kpc for models 
B3, C2, and C3. This demonstrates that more distant encounters excite spiral features in 
regions with larger R. The arms become progressively weaker toward the disk center since 
the ratio of the tidal perturbations to the background gravity is proportional roughly to 
R 2 for small radii. They eventually attain vanishingly small amplitudes inside R ~ 4 kpc 
corresponding to the ILR through which stellar spiral waves cannot propagateQ Although 
the tidal perturbations are strong in extended disks, on the other hand, F still decreases 
with increasing R (> -R max ). This is because the amount of mass available to construct 
spiral arms in the background stellar disk declines very rapidly with radius. Figure [11] also 
shows that self-gravitating spiral arms in model A2 are stronger by about a factor of 1.5 and 
located relatively closer to the center than non-self-gravitating arms in model A2*. 



To study how rapidly spiral arms grow and how long they survive, we plot in Figure 
the temporal variations of F averaged over a range of radii where the arms are strongest 
in each model. Obviously, the arms grow earlier and more rapidly for models with stronger 
tidal perturbations. For instance, it takes only ~ 0.1 — 0.3 Gyrs for the strong-encounter 
models Al and A2 to achieve the peak strength, while more than 1 Gyrs are required for 
the weak encounter models. Figure [13] plots the peak value F max of the arm strength as a 
function of the tidal strength 5", showing roughly F max = O.TQS' ' 83 . 

Since the formation of tidal spiral arms in a disk involves the gathering of particles from 
different radii, the velocity dispersions increase as the arms grow. In addition, gravitational 
scatterings of stellar particles off the arms become efficient to heat the disk once the arms 
acqui re considerable amplitudes , counterbalancing the arm-amplifying e ffect of self-gravity 
(e.g., ISellwood &: Carlbergj|l984l ; iBinney fc Tremaind 119871 ; lBinneyll200ll ). In all the models 



we have considered, the arms never become fully self-gravitating. They stop growing and 
decay as the enhanced velocity dispersions make the once well-organized epicycle orbits 
kinematically less coherent. Figure [T2l shows that for the self-gravitating models, F decreases 



7 It is unclear whether the absence of spiral arms at R < 4 kpc in our models is mainly due to the ILR 
barrier or just because the tidal perturbations are too weak to excite density waves there. We have run a 
model simulation (not listed in Tabic [1]) corresponding to model A2 but without the b ulge (hence no ILR) 
only to find that the inner disk is contaminated by the formation of a central bar (e.g., lNoguchilll987l ). 
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after the peak almost exponentially in a time scale of ~ 1 Gyrs, whereas spiral arms in the 
non-self-gravitating model A2* decay much more slowly since they do not experience secular 
disk heating. Strong encounter models possess spiral arms with F > 5% for ~ 1 Gyrs, 
corresponding to four disk rotations at R = 10 kpc, while spiral arms in the weak-encounter 
model C3 have F £ 3% throughout the entire evolution. Small- amplitude fluctuations of F 
at t £j 1 are caused by the interactions with the particles once pertaining to the bridge and 
tail. 



4.3. Pattern Speed 

Finally, we discuss the pattern speed of the tidal arms formed in our simulations. To 
measure the pattern speed at a given radius, we define the normalized cross-correlation of 
the perturbed surface densities at two fixed times separated by At as 

C ( R M = ^T^2 f 8E(R,<f>,t)SE(R,<f> + 6,t + At)d<f>. (9) 
^o{R) 2 Jo 

For a sufficiently small value of At, the instantaneous arm pattern speed at a given radius 
is then determined by ft p (R,t) = 9 mSuX /At, where # max denotes the phase angle at which 
C(R,6,t) is maximized. We take At = 0.1 in calculating Q p from the numerical data. 

Figure [141 plots as contours the amplitudes of C(R,9,t) on the radius (R)— frequency 
(6/ At) domain for some selected time epochs of models A2 and A2*. The solid and dashed 
lines draw the radial variations of Q and fl± k/2, respectively, from the initial disk rotation. 
At t = 0.1, the spiral arms in both models are relatively weak and the cross-correlation is 
dominated by the extended-disk features, especially by the tidal bridge. The bridge rotates 
almost rigidly at a fixed pattern speed (~ 9.5 km s _1 kpc^ 1 ), corresponding to the angular 
frequency of the perturber at the pericenter. This evidences the phase locking of particle 
orbits in the bridge explained in §3.11 The tail at the opposite side of the perturber becomes 
strong at about t = 0.2, significantly contributing to C(R, 9, t) at R £ 17 kpc. Interestingly, 
the instantaneous pattern speed of the tail is similar to that of the bridge at this time. As 
time evolves further, the extended tidal structures become weaker since the perturber moves 
farther away, while the spiral arms become more pronounced in the distribution of C(R, 8, t). 

When the arms are quite strong (t ~ 0.2 — 0.6) in model A2, their patten speed decreases 
with radius, indicating that they are not a "pattern" in a strict sense. This is the reason 
why the pitch angle of the arms decreases with time. Since the axisymmetric background 
state of the stellar disk as well as the shape and pitch angle of the arms are already known, 
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one can calculate the theoretical pattern speed predicted from the WKB dispersion relation 



2nGE\k R \ , , . 

v = 1 5 Fv\x), (10) 

K 



for tightly- wound density waves (ILin. Yuan, fc Shulll969l ). The dotted line in each of the 
left panels of Figure [TH shows Q p obtained from equation (TTDl) . which traces the loci of 
maximum C(R,6,t) fairly well. Note that equation fllOp would simply yield v ~ — 1 or 
fi p = Q — k/2 without self-gravity, in excellent agreement with the pattern speed of spiral 
arms in model A2* for t £ 0.3. Although the presence of self-gravity tends to enhance 
the arm pattern speed, our numerical results suggest that its effect is quite small; for all 
the models considered, Q p is below ~ 20 km s _1 kpc -1 even when the arms reach the peak 
strength, and it comes very close to the Q — k/2 curve at t 0.6. This implies that the spiral 
arms at least at late time are kinematic spiral waves in which the large velocity dispersions 
of particles as well as the kinematic winding of the arms make self-gravity unimportant. 



5. Summary & Discussion 
5.1. Summary 

Galactic spiral shocks and their substructure-forming instabilities in disk galaxies are 
strongly affected by stellar spiral arms that are often triggered by tidal interactions with 
a companion galaxy. To gain an insight on the large-scale star formation occurring in 
the gaseous component and related evolution of disk galaxies, it is crucial to understand 
the physical properties of tidally-induced stellar arms. While the literature abounds with 
studies of tidal interactions of galaxies, most of them concentrate mainly on morphological 
transformation, especially in the extended parts, of disk galaxies. 

In this paper, we have initiated numerical iV-body experiments for tidal encounters to 
quantify the properties of spiral arms that form in the disks and study how their properties 
vary with tidal strength. We also study the nature of the tidal bridge and tail that develop 
in the outer regions. We consider a simple galaxy model consisting of a rigid halo/bulge and 
a razor-thin stellar disk with Toomre stability parameter of Q ~ 2. A perturbing companion 
galaxy is treated as a point-mass potential moving on a prescribed, prograde, parabolic orbit 
in the same plane as the galactic disk. By varying the mass and pericenter distance of the 
perturber, we explore tidal interactions with strength in the range of 0.03 ^ S £ 0.3, where 
S is the dimensionless momentum applied by the perturber to stars at outer disks (see eq. 
0). 

Our main results are summarized as follows. 
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1. The tidal bridge forms at the near side to the perturber as particles in outer disks 
are pulled out by strong tidal perturbations. Some particles with velocities exceeding the 
escape velocity become unbound, and either are captured by the perturber or escape from 
the system, but these are less than 7% of the total for S ^ 0.3. On the other hand, bound 
particles with low velocities in the bridge execute coherent forced oscillations in such a way 
that the maximum radial velocities Vr and the maximum gradient of the azimuthal velocities 
|9t^/<90| are always attained in the direction toward the perturber. This phase locking of 
the perturbed particle orbits allows the bridge to construct a transient pattern that corotates 
with the perturber as long as the perturber remains close to the pericenter (t £ 0.3). The 
phase locking is also a cause of the sharp leading edge of the bridge, where particles begin 
to fall radially inward during their forced oscillations. 

2. Only strong tidal encounters with S > 0.07 produce a recognizable tail (or counter- 
stream) at the far side of the disk. The tail develops as strongly-perturbed, near-side particles 
overtake mildly-perturbed, far-side particles. When the tail achieves a peak strength, it is 
very narrow and in a roughly logarithmic shape. For 0.07 £j S £ 0.3 we have considered, the 
formation epoch t ta ii, pitch angle i ta ii, and the surface density £ ta ii of the tail depend on the 
tidal strength parameter S as t tail = 0.07S" - 54 , tanz taU = 0.75S - 22 , and E tail /E d = 79S 72 
at R = 20 kpc. Comprising of particles collected from a wide range of radii in the unper- 
turbed disk, the tail is a material arm and has large velocity dispersions, so that it widens 
and weakens with time. 

3. Even though the boost of epicycle amplitudes due to tidal perturbations is quite 
small in regions with R ^ 15 kpc, the perturbed particle orbits are kinematically organized 
to generate two-armed global spiral arms there. With Q ~ 2 in the unperturbed disk, the self- 
gravity of stars does not play a dominant role in growing the spiral modes, although it appears 
to enhance the amplitudes considerably when the arms are nonlinear. The spiral arms are 
approximately logarithmic in shape and subject to kinematic winding. For the parameters 
we have explored, the pitch angle of the spiral arms is in the range of tanz ~ 0.3 — 0.4 when 
the arms attain peak amplitudes and then decreases as tanz oc t-°- 5 ~-°- 6 ; with a smaller 
decay rate corresponding to stronger arms. 

4. Stronger encounter models tend to develop stronger spiral arms earlier and more 
toward the galaxy center, resulting in the arms at R ~ 5 — 10 kpc, ~ 0.1 — 0.3 Gyr after 
the pericenter passage for models with 5* > 0.13. Arms are absent inward of R = 4 kpc 
corresponding to the inner Lindblad resonance. In terms of the parameter F (eq. [8]) that 
measures the perturbed radial force due to the spiral arms relative to the mean axisymmetric 
gravity, the maximum strength of the spiral arms behaves as -F max = 0.79S 10 ' 83 . Because of 
large velocity dispersions associated with the particle gathering and secular heating, the 
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arms never become fully self-gravitating and decay after the peak almost exponentially in a 
time scale of ~ 1 Gyr. 

5. Analyses using the normalized cross- correlation of the perturbed densities reveal that 
the arm pattern speed fl p is not constant in both radius and time, indicating that spiral arms 
that form in our models are not exactly a pattern. In fact, Q p decreases with radius, causing 
the pitch angle to decrease with time. Self-gravity tends to increase Q p , but only below 
~ 20 km s _1 kpc -1 even when the arms are strongest. Self-gravity becomes unimportant as 
the arms decay, resulting in Q p « Q — k/2 at late time. 



5.2. Discussion 



We have seen in this paper that spiral arms produced by tidal encounters are approxi- 
mately logarit hmic in shape, s i milarly to obseryed spiral arms in many g rand-design spiral 
galaxies (e.g., lKennicuttlll98ll ; lElmegreen et al.lll989l ; IShetty et al.l 120071 ). The occurrence 
of the logarithmic arms in our models can be understood as follows. As mentioned above, 
the arms are kinematic density waves modified by self-gravity. Ignoring the effect of self- 
gravity and assuming that the phases of the waves are aligned along (p = (f> p = at t = 0, 
corresponding to the impulsive tidal perturbations applied at the pericenter, the pitch angle 
of kinematic density waves with m = 2 is given by equation (J5J). If the right-hand side of 
equation (jSJ) is independent of R, the arms have a perfect logarithmic shape. It turns out 
that the galaxy rotation curve we adopt (Fig. [T5|) has an approximately constant value of 
d(Q — k/2)/ din R ~ 3.5 ± 0.5 km s _1 kpc -1 over the distance from the ILR radius out to 
the edge of the disk. This results in |Atanz|/ tani ~ 0.15 over a range of radii where spiral 
arms are strong, indicating that the variation of the pitch angle along the arms is in fact 
very small. The presence of self-gravity as well as epicycle motions of particles are likely to 
further smooth out the local variation of tani. 

Our numerical results show that self-gravity is unable to keep the arm pitch angles 
fixed over time. A larger rate of shear in the rotation cu rve implies a smaller arm pitch 
angle for kinematic arms. Indeed, ISeigar et al.l (120051 . 120061 ) reported a well-defined negative 
correlation between the arm pitch and the shear rate for a sample of (not necessarily tidally- 
driven) spiral galaxies, suggesting that spiral arms in real galaxies are unlikely to be fully 
self-gravitating. 

While we adopt highly simplified models for both the disk galaxy and the orbital param- 
eters of tidal interactions, it is still interesting to compare the arm properties found in our 
simulations with those of observed spiral arms. In the case of the M51/NGC5195 system, 
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the mass ratio of t he target galaxy to the com panion is estimated to be ~ 0.3 — 0.55 (e.g., 
Smith et al.l Il990l ; ISalo &: Laurikainenl l2000al ). The encounter models that well reproduce 
the kinematics and morphol ogies of the M51 system fav or an inclined orbit with the percen- 
ter distance of 20 — 30 kpc (ISalo &: Laurikainenl l2000al ). Since the thin disk approximation 
and non-inclined orbits taken in our models tend to produce stronger tidal arms than in the 
thick-disk, inclined-orbit counterparts, models Al and A2 can perhaps be best compared 
with the M51/NGC5195 system. K-band observat ions indicate that the radially-averaged 
spiral arm stre ngth F is around 2 % for M51 (e.g. . 



2000bl ; see also iRix fc Riekel Il993l : iRix fc Zaritskv 
~ 17 



Scoville et al. 2001; Salo fc Laurikainen 



1993), which is not much different from 
22% found for models Al and A2 at t ~ 0.1 - 0.3 ( Fi gs. [TT1 and [T2[). Th e arms in 
M51 are logarithmic spirals with a pitch angle of tanz ~ 0.39 (jShetty et al.l 120071 ). which is 
again close to the arm pitch angle in model A2 at t ~ 0.2 — 0.3. 



Among the properties of spir al arms, the most i ntrigu ing is the pattern speed that is not 
well constrained by observations. lElmegreen et al.l (119891 ) identified 4 : 1 resona nce features 
in the arms of M51 to find Q p ~ 40 km s _1 kpc -1 , while IZimmer et al.l (120041 ) determined 
Vl p = 38 ± 7 km s" 1 kpc -1 using the Tremaine- Weinber g method. By running col lisional 
models for cloud dynamics under a given spiral potential, iGarcia-Burillo et al.l (119931 ) found 
Q p ~ 27 km s _1 kpc -1 for the best fit to the observed morphologies of the CO arms in M51. 
All of these works were based on the premise that the arm pattern speed is a constant with 
radius. However, our numerical results show that the pattern speed of tidal arms depends 
on the radius. In the case of model A2, fl p is a decreasing function of radius, varying when 
the arms are strongest from ~ 20 km s _1 kpc^ 1 at the ILR to ~ 10 km s _1 kpc -1 at the 
ou ter parts, and at late r time c onverging to the Q — k/2 curve. A similar trend was obtained 
by ISalo fc Laurikainenl (l2000bl ) who ran more realistic encounter models (with a star-only 
disk) for the M51 system and found that fl p is close to the Q — k/2 curve for a range of radii 
where the spiral arms are strong. Although much remains uncertain regarding the effects 
of the cold gaseous component and rotation curve, these results suggest that tidally- driven 
arms may have a pattern speed that varies with radius in real spiral galaxies. 

An age distribution of star c lusters in M5 1 shows a narrow peak at 4 — 10 Myrs and 
a broad peak at 100 — 400 Myrs (ILee et al.l 120051 ). indicating active star formation at these 
epoches. This enhanced star formation is most likely due to strong spiral arms induced by the 
tidal interactions with the companion NGC 5195. Since it takes about ~ 100 — 200 Myrs for 
the spiral arms in our models Al and A2 to attain a substantial amplitude, say F = 10%, 
after the perturber passes the pericenter, this implies that the closes t passages of NGC 
5195 m ight have occurred ~ 100 — 200 Myr and ~ 200 — 600 Myrs ago. ISalo fc Laurikainen 
(j2000al ) proposed two encounter models for the M51 system: a near-parabolic, single-passage 
orbit occurred 400-500 Myrs ago and a bound double-passage orbit having taken place 400- 
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500Myrs and 50-100 Myrs ago. Considering the delay between the pericenter passage and 
the development of strong arms, the cluster age distribution appears to be more consistent 
with the double-passage scenario, although it is uncertain what effects the second passage 
will make on the pre-existing arms generated at the first passage. 



It is well known from the seminal paper of iToomre fe Toomrei (119721 ) that tidal inter- 
actions distort the outer parts of a galactic disk and create a tidal bridge extending toward 
the perturber as well as a narrow tail at the opposite side. They noted a fraction of the disk 
material is stripped and transferred through the bridge to the perturber. In this work, we 
further show that the bridge is in fact a transient pattern constructed by bound particles 
whose orbits are strongly locked to the perturber. As these particles follow galactic rotation, 
they are pulled out toward the perturber and then move radially inward at the leading edge, 
making the bridge rather sharp. By mapping the final to initial particle positions under 
an impulse approximation, iDonner et al.l (119911 ) showed that the sharp boundary of a tidal 
bridge corresponds to the loci (caustics) of zero Jacobian of the mapping where the orbits 
of neighboring particles come very close together. Indeed, Figure H] shows that the leading 
edge has a large velocity dispersion, consistent with the Liouville theorem that dictates the 
conservation of the particle density in the phase space. 

Unlike a bridge, a tail at the opposite side is a material arm resulting from the overlap- 
ping of near-side particles with far-side particles in the extended parts of the disk. Conse- 
quently, t he tail forms l a ter th an t he bridge by abou t a hal f orbital time, consistent with the 
results of IDonner et al.l (ll99ll ) and iByrd fe Howard! (119921 ). Our experiments show that the 
formation time and pitch angle of a tail are well correlated with the tidal strength parameter 
S. While we employed simple models for tidal interactions and limited our simulations to 
the cases with S < 0.3, our results appear to be applicable to models wit h quite s t rong tidal 
perturbations as well. In simulations of merger encounters, for example, iBarned (119921 ) ran 
self-consistent models consisting of a live halo/bulge and a disk with both stars and gas. One 
of his models considered interactions between equal-mass disk galaxies, in which one disk 
passes directly through t he other with the pericenter distance R per i/R g = 0.5, corresponding 
to S = 1.48. Figure 3 of IBarned (119921 ) shows that the tail in this model becomes strongest 
at t ~ 1.25, corresponding in our units to t ta ii ~ 0.053 after the pericenter passage, and has 
a logarithmic shape with tani tai i « 0.83, which are remarkably similar to the extrapolation 
of our results in §3.21 that yiel d tta.ii ~ 0.057 and tanitaij ~ 0.81. Through a comprehensive 
survey of the parameter space, IToomre fe Toomrei (119721 ) found that tail shape is insensitive 
to the orbital eccentricity e for 0.6 < e < 1 as long as the inclination of the orbit is not so 
large (see also IBarned Il998l ) . which is also consistent with our result that tanz ta ii is weakly 
dependent on S. 
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Numerical studies on tidal encounters often report t he formation of double arm structure 
at th e opposite side to the perturber (e.g.. ISundinlll989l ; lElmegreen et al.lll991t iDonner et al. 
199ll ). Our simulations also exhibit such double features (see, e.g., t = 0.3 frame of Fig. [T]) 



which come out as the tidal tail decouples from the spiral arms tha t, because of the smooth 
alignment with the former, are not readily discernible at t = 0.2. lElmegreen et al.l (jl99fl ) 
found that the lagging arm forms by gathering particles streaming away from the near side 
and soon merges with the leading arm. This might be a consequence of the zero velocity 
dispersion in their unperturbed disk since the ratio of the velocity impulse due to tidal torque 
to the in itial random veloc i ty is t oo large to set up well-defined spiral arms in the disks of their 
models. Elmegreen et al. fjl99lh also found that a prograde, in-plane encounter produces a 
"ocular" galaxy with oval-shaped, sharp boundaries, provided S > 0.019. A similar structure 
can be seen in the t = 0.2 panel of Figure dj although the boundaries in our models are less 
sharp since, as they noted, the formation of ocular shape requires the injected energy from 
the perturber to be much larger than the kinetic energy in random particle motions. 
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A. Galaxy Model 

In this Appendix we describe the model galaxy we use for tidal encounter experiments. 
The galaxy consists a rigid halo/bulge and a live stellar disk. For a fixed spherical halo, we 
adopt a truncated logarithmic potential 

(T, ( \ — / \ v o ( r c + r2 ) + constant forr < r tr 
h{r) ~ { -GM h (r tI )/r forr > r tr ( } 

where r is the three-dimensional distance from the halo center, r c is the halo core radius, 
rtr is the truncation radius, and vq is the cons tant rotation vel ocity the disk would have 



at large r if the halo were not truncated (e.g., iLee et al.J Il999l ). The corresponding halo 
mass distribution is M h {r) = v$r 3 /[G{r 2 + r 2 )} for r < r tr and M h (r) = M h (r tl ) for r > r tr . 
The constant in equation (I All) should equal — v^rlJir 2 + r tr) ~~ h v o l°s( r c + r tr) to make 
the potential continuous at r = r tr . For the simulations presented in this paper, we take 
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r c = 7.5 kpc, r tT = 25 kpc, v = 220 km s 1 , corresponding to M/ l (r tr ) 
spherical bulge is modeled by a Plummer potential 



2.58 x 10 11 M . A 



GM h 



\Jr 2 + a 2 



(A2) 



with the scale radius a = 0.23 kpc and the total bulge mass M& = 1.0 x 10 M Q . 

Although stars in real galactic disks are distributed with a finite vertical thickness 



for example, amou nting; to ~ 330 pc in the solar neighborhood (e.g., iChen et al.l 12001 
Karaali et al.l 120041 ). we impose an infinitesimally-thin stellar disk by setting the vertical co- 
ordinates and velocities equal to zero throughout the simulations. For the radial distribution 
of stellar surface density, we adopt an exponential form 



X i (R) = E exp(-R/R d ), 



(A3) 



where R is the galactocentric radius in the disk, R d is the disk scale length, and S is the 
surface density at the galaxy center. The total disk mass is M d = 2tiT, q R 2 ,. The gravitational 
potential of the disk is given by 



$ d (R) = ~{GM d /R d )R I Q {R)K X {R) - h(R)K (R) 



(A4) 



where I n and K n represent modifi ed Bessel functions of the first and second kinds, respec- 
tively, and R = R/2R d (see e.g., iBinney &: Tremaind 119871 ). We take R d = 3.4 kpc and 
S = 711 M pc~ 2 , corresponding to M d = 5.2 x 10 10 M . 



To obtain the equilibrium velocity distribution of disk particles u nder the to t al gra vita- 
tional potent i al $ to t = $h + ^6 + we follow a method suggested by iHernquistl (119931 ) and 
Quinn et al.l (I1993I ). We first assume that the radial and azimuthal components, vr and v^, 
of particle velocities obey initially the Schwarzschild distribution function 



Iixoro^ 



exp 



Vr (V,/, - V^f 



2a 2 R 



2ai 



(A5) 



wher e <j R and cr^ are the radial and azimuthal velocity dispersions, respectively (e.g.. lToomre 
19641 ). The mean azimuthal streaming velocity differs from the circular velocity v c deter- 
mined solely from the total gravitational potential as v 2 (R) = —d<$> tot /dhiR. In the local 
approximation in which S^, <Jr, and are assumed to vary slowly with R, one can show 
that or and a a, are related to each other through 



al/al = n 2 /An 2 



(A6) 
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where Q = f c /i? is the local rotational a ngular velocity and k 2 = 4ft 2 + dQ 2 /d In i? is the 



square of the local epicycle frequency (e.g.. lBinney fc Tremaindll987l ). Then, the usual Jeans 



equation in the radial direction for an equilibrium disk leads to 



en 



1 - 



4ft 2 



2* 

Rd 



(A7) 



f lBarneslll992l : lHernquistlll993f ). 



Finally, we express the radial velocity dispersion or in terms of the Toomre stability 
parameter 

« = ife (A8) 

which determines local gravitational stability of a razor-thin disk to axisymmetric perturba- 
tions. We adopt a fixed value of Q = 2 everywhere initially. This value of Q corresponds 
roughly to solar neighbor hood conditions with k = 36 km s~ l kpc" 1 . Or = 30 km s~ x 
( iBinney fc Tremaind 119871 ). and = 35 M Q pc -2 ( jKuijken fc Gilmoreill989l ). and is large 



enough to make swing amplification of non-axisymmetric disturbances inefficient. This pre- 
cludes the possibility of spiral arms driven sp ontaneously by the stellar self-gravity (e.g., 
Sellwood & CarlberelE983lBertin et alill989bh . 



Figure [15] plots the circular velocity v c (R) and the mean rotational velocity v^R) of 
our model galaxy as solid and dotted lines, respectively. Also shown as dashed lines are 
the separate contributions to v c from halo, bulge, and disk, which have a mass ratio of 
Mh : Mb : = 0.81 : 0.03 : 0.16 inside R = 25 kpc. It is apparent that is usually 
smaller than v c , indicating that stars, on average, lag behind a circular orbit at the same 
galactocentric radius, a phenomenon known as asymmetric drift. 



B. Initial Disk Setup 

We initialize the exponential stellar disk (eq. |A3j ) by distributing iV=514,000 equal- 
mass particles and place it under the combined halo and bulge potential (eqs. |Alj and 
|A2| ). Strictly speaking, the model disk constructed in this way is not in perfect equilibrium 
because equations (1A5I) and (1A7j) hold true only in a local sense, that is, only when the 
gravi tational potent ial and the stellar velocity dispersions do not vary much with radius 



(e.g. 



Sellwoodlll985l ). In addition, when the disk is allowed to evolve, any non-axisymmetric 
modes that grow may interact with particles, feeding them with random kinetic energy. 
Two-body interactions of particles tending to heat the disk are not completely negligible, 
either. All of these may cause the disk structure to deviate considerably from the desired 
one even before undergoing tidal encounters. 
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To obtain a disk configuration representing a dynamically well-relaxed, global equilib- 
rium, we evolve our model galaxy in isolation for two Gyrs. Figure [16] displays snapshots of 
particle distributions from the isolated disk evolution. The disk is rotating in the counter- 
clockwise direction and time is expressed in units of Gyr. Other than weak non-axisymmetric, 
trailing structures seen at its outskirts, the disk does not suffer from dramatic morphological 
changes. This implies that the disk is globally stable, a consequence of the fact that, when 
Q ~ 2, the growth of perturbations by swing amplification and other instabilities is quite 
mild (jToomrd Il98ll ; ISellwoodl Il989l ) . No ad ditional perturb ation from the rigid halo and 
bulge also helps to keep the disk featureless (IHernquistl Il993l ) . 

Figure [T7| shows the radial distributions of various physical quantities averaged over the 
azimuthal direction at t — 0, 1, and 2 Gyrs. While v$ and gr change promptly (within less 
than 0.1 Gyr) from the initial profiles, remains almost unchanged. The changes in gr 
and Q are largest at R ^ 5 kpc where the circular velocity (hence the total gravitational 



potential' 


varies 


Sellwood 


1985). 



> 



10 kpc from the initial values are likely 
caused by mild swing amplification. Except the slight variations of gr near the center, the 
changes of the disk properties between 1 and 2 Gyrs are practically negligible, indicating 
that at late time the disk is in a sufficiently well-relaxed, new equilibrium. 



C. Impulse Approximation 

In the absence of tidal perturbations, the motions of individual disk particles are in 
general a superposition of the radial oscillations with epicycle frequency k around their 
guiding centers and the circular oscillations of the guiding centers about the disk center. The 
dispersion 5R in the epicycle amplitudes is related to the radial velocity dispersion through 
5R = gr/k. Tidal perturbations are able to enhance the epicycle amplitudes for particles 
whose orbital periods are not so small compared with the duration of a tidal encounter. 

Using an impulse approximation, one can estimate 5R of disk particles subject to tidal 
perturbations. Let us assume that the tidal forcing is applied impulsively near the pericenter 
during the time interval of R per \/v p . Then, the increment Avr in the radial velocities of 
particles at radius R is given approximately by 

i A i 2GM p Rq 

\ AV R\ = ~ D2 » ( C1 ) 

wher e v p = [2G(M^ + M^lMpsr i] 1S the orbital velocity of the perturber at the pericenter 



e.g.. lBinney fc Tremaind 119871 ). Assuming that the kinetic energy associated with \ Avr\ is 
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absorbed into the epicycle motions, one obtains 

5R = (a| + AvrY^/k (C2) 

as a measure of the mean radial excursion of disk particles under the influence of tidal 
perturbations. Figure [3] plots as thin curves 5R from equations (I CI I) and (1C2I) with differing 
M p corresponding to models A2, B2, and C2, while the thick curve draws (Tr/k. 
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Table 1. Summary of model parameters and simulation results 



Model a 


Mp/M g 


-Rperi(kpc) 


S b 


Mcap/M d c (%) 


M csc /M d d (%) 


*t ail 


tan i tai i 


£tail/£20 e 


^max 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Al 


0.4 


25 


0.250 


4.99 


1.97 


0.14 


0.546 


29.2 


0.22 


A2 


0.4 


35 


0.151 


1.99 


0.08 


0.18 


0.511 


22.0 


0.18 


A2* 


0.4 


35 


0.151 


1.84 


0.08 


0.19 


0.707 


11.6 


0.11 


A3 


0.4 


45 


0.103 


0.46 


0.006 


0.25 


0.499 


15.3 


0.09 


Bl 


0.2 


25 


0.135 


1.95 


1.33 


0.18 


0.457 


17.4 


0.17 


B2 


0.2 


35 


0.081 


0.79 


0.03 


0.24 


0.411 


13.4 


0.10 


B3 


0.2 


45 


0.056 


0.07 


0.002 








0.06 


CI 


0.1 


25 


0.070 


0.73 


0.67 








0.11 


C2 


0.1 


35 


0.043 


0.22 


0.01 








0.06 


C3 


0.1 


45 


0.029 












0.04 



a Modcl A2* is identical to model A2 except that the former neglects the self-gravity of the perturbed density. 

b S = (M p /M g )(R g /R peii ) 3 (AT/T) = 0.738(M p /M s )(iJ 9 /R p ) 3 / 2 [l + (Mp/Mg)]- 1 / 2 is the dimcnsionlcss tidal strength parameter. 
°A^cap is the total mass of the captured particles by the perturbing companion. 
d M eS c is the total mass of the non-captured, escaping particles from the whole system. 
e S2o = S d (20 kpc) is the surface density of the initial unperturbed disk at R = 20 kpc. 
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Fig. 1. — Snapshots of the particle distributions in model A2 in the x-y plane. Only 20% of 
the particles are shown to reduce crowding. The elapsed time is shown at the upper right 
corner of each panel. The arrow and the associated number give the direction and distance 
(in kpc) to the perturber that passes through the pericenter (x,y)—(35 kpc, 0) at t = in 
the counterclockwise direction. See text for details. 
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Fig. 2. — Distributions of the perturbed surface density of model A2 in the polar 

coordinates. When t ^ 0.2, is dominated by the extended-disk structures such as bridge 
and tail at R ^ 15 kpc, which become loose and spread widely after t ~ 0.3. The spiral arms 
at R £ 15 kpc appear straight in the — lni? plane, with a slope becoming progressively 
smaller with time. 
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Fig. 3. — Dispersions in the radial departures of particles from the initial locations over the 
course of tidal interactions. Various symbols indicate the numerical results of models A2, B2, 
and C2. The thin solid, dotted, dashed curves draw the analytic estimates (eq. |C2j ) based 
on the impulse approximation. The thick solid line corresponds to the unperturbed disk in 
which 5R is purely due to the epicycle orbits associated with the initial velocity dispersions. 
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Fig. 4. — Azimuthal variations of the radial (vr) and azimuthal (v<p) velocities of the particles 
at R = 20 kpc in model A2 at early epochs of tidal interactions. The azimuthal phase in the 
abscissa is repeated for clarity. In each panel, vertical dotted lines mark the phase angles of 
the perturber. 
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Fig. 5. — Spatial distributions of selected particles in model A2 at some time epochs when 
the perturber is close to the pericenter. Blue, green, yellow, orange, and red colors represent 
the particles originally located in annuli with i?o = 6 — 7, 10 — 11, 14 — 15, 18 — 20, and 
22 — 24 kpc, respectively. The black circle in each panel has a radius of 20 kpc and the arrow 
indicates the direction to the perturber. The black and cyan dots denote the groups of the 
near-side and far-side particles, respectively, in the Ro = 22 — 24 kpc ring, which merge at 
the far side to form a tail at t — 0.2. 
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Fig. 6. — Dependences on the the tidal strength parameter S of the formation epoch t ta ii, 
the pitch angle ? tail at t — t ta ih an d the surface density S tai i at R = 20 kpc and t = t ta n of 
tidal tails in various models. The dotted line in each panel gives the best fit to the numerical 
results: t tai i = 0.07S' - 54 , tanitaa = 0.75S - 22 , and S tai i/S 2 o = 79.0S - 72 , where S 20 denotes 
the surface density of the initial disk at R = 20 kpc. 
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Fig. 7. — Close-up views of the distributions of stellar surface density in model A2, with 
the color bar labeling (E/10 4 M pc -2 ) 1 / 2 . The two-armed spiral density waves excited by 
the tidal forcing achieve the maximum strength at t ~ 0.3, and then gradually weaken. 



-40- 




Fig. 8. — Evolution of the Fourier amplitudes with wavenumber p, defined by equation (jlj), of 
the m = 2 logarithmic spirals in model A2. The modal growth is due to swing amplification 
at early time (t £ 0.04), which becomes soon dominated by the kinematic overlapping of the 
perturbed epicycle orbits. 
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Fig. 9. — Temporal changes of the pitch angle of the spiral arms located in the R — 5 — 10, 
8.0-8.5, 8-13, and 11-16 kpc regions for model A2, A2*, B2, and C2, respectively. In the 
self-gravitating models A2, B2, and C2, the pitch angle deceases as tan 2 oc t~ a5 a6 , with 
weaker arms decaying slightly more rapidly, whereas tani oc t' 1 for the non-self-gravitating 
model A2*. 
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Fig. 10. — Azimuthal distributions of the perturbed density <5£ (top) as black curves and 
the radial velocity vr (middle) and perturbed circular velocity dv^ = v $ — (bottom) as 
dots in the disk at R — 10 kpc when t = 0.4 for model A2 (left) and when t = 0.5 for model 
C2 (right). Note that the vertical scales are different from the left and right panels. In the 
top panels, the red lines give the m = 2 Fourier modes 5E m=2 of the perturbed density. In 
the middle and bottom panels, the blue curves give the average values of vr and Sv^, while 
the red curves plot the predictions of the linear density wave theory corresponding to <5£ m=2 . 
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Fig. 11. — Arm strength F averaged over the time interval At = 0.4 centered at the time of 
the peak strength as a function of radius. The induced spiral arms for stronger encounter 
models Al, A2, Bl, and CI peak at -R max ~ 5 — 10 kpc, while weaker encounter models B3, 
C2, and C3 produce spiral arms at -R max ~ 11 — 16 kpc. 




Fig. 12. — Time evolution of the arm strength F averaged over 5 kpc R 10 kpc for 
models Al, A2, Bl, and CI, over 8 kpc <, R 13 kpc for models A2*, A3 and B2, and over 
11 kpc R £ 16 kpc for models B3, C2, and C3. In each self-gravitating model, it takes 
about one or two rotational periods at R = -R max for the arms to reach the maximum value. 
After the peak, F decays as ~ exp(— 1/\ Gyr) due largely to large dispersions in the particle 
velocities. Fluctuations of F at t ^ 1 arise as the particles once pertaining to the bridge and 
tail move in and out the arms. Without disk heating, the decay of the spiral arms in the 
non-self-gravitating model A2* is quite slow. 
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Fig. 13. — Dependence of the peak arm strength -F max on S. The dotted line F max = 0.955 10 ' 
is the best fit to our numerical results. 
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Fig. 14. — Contours of the cross correlation of the normalized surface density in the radius 
- frequency domain for models A2 (left) and A2* (right). Smooth curves draw Q (solid) 
and Q ± k/2 (dashed) from the initial disk rotation. The dotted line in each of the left 
panels plots the theoretical patten speed calculated from the linear dispersion relation for 
the background parameters equal to the azimuthally-averaged disk values obtained from the 
simulation. Nearly constant tt p ~ 9.5 km s" 1 kpc -1 at R ^ 17 kpc for t £ 0.2 traces the 
tidal bridge and tail, which instantaneously corotate with the perturber. The loci of the 
maximum cross correlation for the arms match well with the Vt — k/2 curve in model A2 for 
t £ 0.6 and in model A2* for t £ 0.3. 
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Fig. 15. — Circular speed v c (R) (solid) and mean streaming velocity v^(R) (dotted) of our 
model galaxy as functions of the galactocentric radius R. Contributions to v c from disk, 
bulge, and halo are plotted as dashed lines. 
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Fig. 16. — Evolution of our model disk in isolation. Time is expressed in units of 10 9 yr. Only 
10 4 particles are plotted to reduce crowding. The disk is rotating in the counterclockwise 
direction. No notable change in appearance is found, indicating that the disk is globally 
stable. 
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Fig. 17.— Radial distributions of (a) the disk surface density E^, (b) mean streaming 
velocity v^, (c) radial velocity dispersion ctr and (d) Toomre stability parameter Q from an 
isolated disk evolution at t — (dotted), 1 Gyr (solid), and 2 Gyr (dashed). 



